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INTRODUCTION 


This  final  technical  report  describes  work  carried  out  on  the  Bio-Inspired  Concepts  grant 
titled  “Studies  of  Biological  Response  to  External  Electric  Fields  for  Cellular 
Manipulation  and  Diagnostics  -  Modeling  and  Experimentation.”  Significant  progress 
(both  in  experimental  activities  and  in  theoretical  modeling)  has  been  made  at  Old 
Dominion  University,  in  conjunction  with  Eastern  Virginia  Medical  School,  and  is 
described  in  detail  in  this  report. 


LIST  OF  OBJECTIVES  AND  ACCOMPLISHMENTS 

Objective  1:  Modeling  for  Cellular  Response  to  External  Electric  Fields 

Appropriate  bio-models  for  the  analysis  of  cellular  response  to  high-intensity,  short- 
duration,  external  electric  fields  were  developed.  The  focus  was  on  the  three  distinct 
tasks:  (a)  Development  of  numerical  codes  for  three-dimensional,  time-dependent 
solution  of  distributed  potentials  and  current  flows  in  a  cell.  These  models  can  also  be 
applied  to  the  evaluation  of  current-voltage  characteristics  of  membrane  pores.  Polarity 
dependencies  and  asymmetric  behavior  resulting  for  internal  fixed-charge  distributions 
and  non-uniformities  in  pore  cross-sectional  area  can  be  probed,  (b)  Use  of  the  distributed 
electrical  model  as  the  input  to  predict  cell  irreversibility  and  its  pulse-width  dependence 
based  on  a  simple  two-valley  energy  model,  and  (c)  Development  of  Molecular 
Dynamics  simulation  codes  for  microscopic  analysis  at  the  atomic  level  of  cell 
membranes.  The  dynamics  of  pore  formation,  the  anode-side  preference,  and  mechanism 
of  phosphotidylserine  (PS)  externalization  was  examined 


Accomplishments  Under  Objective  1 

(A)  Development  of  a  Macroscopic.  Distributed  Cell  Model; 

All  of  the  three  above  tasks  were  successfully  carried  out.  Models  were  developed  and 
corresponding  numerical  codes  implemented.  The  results  and  predictions  obtained  were 
published  in  various  scientific  journals  [1-6].  Figure  1  shows  a  schematic  for  the 
distributed,  three-dimensional,  time-dependent  model  developed  for  analysis  of  voltages 
(including  trans-membrane  potentials)  and  current  flows.  A  spherical  simulation  volume 
used  for  the  calculations  is  shown  with  an  axial  electric  field  pulse  encompassing  a  cell 
and  an  internal  organelle.  Pores  have  been  shown  in  this  schematic  as  a  qualitative  guide. 
The  distributed  electrical  representation  of  the  cell  shown  in  Fig.  1,  comprised  of  parallel 
resistive  (R)  and  capacitive  (C)  elements.  Based  on  the  electric  field  and  current 
distributions,  the  power  dissipation  could  be  calculated.  The  maximum  temperature  rise 
within  the  cell  volume  could  also  be  computed.  A  typical  result  for  a  225  kV/cm  voltage 
pulse  having  a  1.5  ns  rise  time,  a  10  ns  ON  time,  and  a  1.5  ns  fall  time  is  shown  on  Fig.  2 
for  a  normal  B-cell.  The  parameters  for  the  B-cell  were  taken  from  the  literature  [7],  A 
maximum  temperature  rise  of  about  1.2  °K  is  predicted.  The  highest  temperature  rise  is 
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predicted  to  occur  just  outside  the  cell,  with  the  cytoplasm  being  a  somewhat  lower 
temperature,  and  the  inner  organelle  having  the  least  heating.  In  actual  experiments,  the 
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Fig.  1.  Schematic  for  the  cellular  calculations,  (a)  The  simulated  region  including  the 
cell  and  an  organelle,  and  (b)  the  distributed  electrical  representation. 
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voltages  (V)  are  typically  on  the  order  of  100  kV/cm.  Since  the  power  dissipation  scales 
as  V2 ,  one  can  expect  a  maximum  heating  of  about  0.3  °K.  This  is  in  line  with  the  notion 
of  an  essentially  non-thermal  biological  effect.  Unfortunately,  it  is  very  difficult  to 
obtain  temperature  measurements  with  temporal  resolution  on  the  nano-second  time 
scales.  Hence,  experimental  corroboration  has  not  been  possible  so  far. 


E  =  225kv/cm,  1.5ns  rise  time  +10ns  on  time+1.5ns  fall  time,  normal  B-cell 
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Fig.  2.  Results  of  the  temperature  rise  in  response  to  an  external  voltage  pulse. 


The  distributed,  macroscopic  electrical  model  of  the  cell  was  used  to  then  used  to 
simulate  poration  at  both  the  outer  (plasma)  and  the  organelle  membranes.  A  continuum 
approach  based  on  the  Smoluchowski  equation  [1,8-12]  was  used  to  model  the  pore 
dynamics.  For  purposes  of  numerical  computations  the  entire  cell  was  broken  up  into 
segments,  and  each  segment  represented  by  a  parallel  RC  combination.  Azimuthal 
symmetry  was  used  to  map  the  three-dimensional  structure  into  the  r  and  cp  co-ordinates 
of  a  spherical  system.  This  method  is  different  from  the  Legendre  polynomial  model 
used  in  the  past  for  electrostatic  cases.  Volume  and  shape  changes  of  cells  were  ignored 
since  the  external  applied  pulse  is  too  short  (~ns)  for  cells  to  deform  much  during  this 
time  interval. 

Time  Dependent  Electrical  Calculations-' 

In  our  simulations,  the  computational  region  is  a  sphere  that  includes  the  cell,  its  sub¬ 
structure,  and  the  surrounding  suspension  medium.  It  is  discretized  in  r  and  #  direction 
as  shown  in  Fig.  3.  For  simplicity,  only  a  quarter  of  the  computational  region  is 
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considered.  For  simplicity,  membranes  are  taken  as  an  integral  unit,  i.e.,  this  sub-region  is 
not  further  discretized.  For  each  element,  the  current  continuity  equation  given  below 
holds: 

V.(J+^)  =  0  ,  (1) 

dt 

where  J  is  the  current  density  and  D  the  electric  displacement  vector.  For  nodes  with 
index  0<i<n  and  0<j<m,  as  shown  in  Fig.  3,  the  above  equation  can  be  rewritten  in  the 
following  form  : 

6  clF  6 

+  )»xA  =Z'«=°  •  (2) 

*=i  dt  *=i 

outer 


membrane  inner 
(~nm)  membrane 


Figure  3:  Schematic  of  one  quarter  of  the  model  used  to  represent  a  cell  for  the 
distributed  electrical  calculations.  The  dotted  box  shows  a  typical  element  with  current 
flows. 


Here  7.5  and  I6  are  currents  along  the  (j)  direction.  Considerations  of  geometric  symmetry 
of  the  computation  region  leads  to  Is  =  I6  =  0  due  to  the  equi-potentials.  Consequently, 
one  obtains  the  following  equation  discretized  in  space  and  time: 
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In  the  above,  V/j  stands  for  the  potential  at  node  (i,  j)  at  time  t.  The  radial  distance 
between  nodes  (i,  j)  and  (n,  m)  is  r.  The  areas  A,  through  A4  are  shown  in  Fig.  3.  Also, 
G|  through  g4,  and  £1  through  £4  are  the  conductivities  and  permittivities,  respectively,  at 
faces  1-4.  Thus,  if  a  node  were  to  be  within  the  environment  medium,  then  all  the  g 
values  would  be  the  same.  However,  for  nodes  on  a  membrane,  G|  would  be  different 
from  02. 

In  order  to  reduce  the  computation  load,  special  boundary  conditions  were  applied, 
and  only  a  quarter  of  the  entire  spherical  computational  region  was  considered  based  on 
the  inherent  symmetry.  Nodes  with  j  =  0  and  j  =  m  had  to  be  treated  carefully.  Only  I,,  12 
and  /4are  nonzero  since,  the  targeted  element  only  has  5  faces  as  face  3  shrinks  to  a  line. 
For  y  =  0  and  0  <  i  <  n,  eqn.  (3)  becomes: 


vLj^vf, 

Ar 

V'  -V' 
v_mj  v_h± 

Ar 

V'  ,  -V' 


(r  +  Ar/2)2  (cos  6  -  cos(0  +  AG!  2))  + 


(r-Ar/2)2  (cos  6  -  cos(0  +  AO  /  2))  + 


(Ja  — — —  rAr  sin(0  +  d9 /  2)  + 


4  rAO 
P  y'+i  _y'+1  y  _i /' 

— - - — - — ](r  + Ar/2)2(-l  -cos(0  + A0/2))  + 


(4) 


At  Ar 

rf+1  1//+1 


Ar 


F  V'  -V' 

— hi - !2hL - i±](r  -  Ar /  2)2 (-1  -  cos(0  +  AO/2))  + 


At  Ar 

7/+1  I//+1 


A  r 


I//+ 1  _i//+l  \rt  _  t/  / 

£±[_hth__jj - hil 1 - — ]rArsin(0  +  AO  12)  =  0 

At  rAO  rAO 


In  the  above,  0  =  n+  jxA0  =  n  and  AO  is  — /m.  For  an  electric  field  applied  applied 

along  the  z  direction  (as  shown  in  Fig.3)  and  potential  at  node  with  /  =  n  is  zero,  i.e, 
Vnj=0.  Nodes  with  index  j  =  0  are  equi-potentials,  so  Vim=Vnj=  0.  Another 
boundary  condition  to  be  considered  is  for  nodes  with  i  =  0.  The  potentials  of  such  nodes 
is  calculated  as:  V0  j  =-E0RcosO  ,  for  j  =  0...m,  where  Eo  is  the  externally  applied 
electrical  field,  and  R  the  radius  of  the  computational  region. 
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Combining  with  the  boundary  conditions  discussed  above,  one  gets  N  equations  for  the  N 
unknown  node  voltages.  These  N  equations  can  easily  be  solved  through  matrix  inversion 
of  the  form:  Ax  =  B.  Here  A  is  a  sparse  coefficient  matrix.  Potential  on  each  node  is 
easily  updated  at  each  time  step. 

A  summary  of  the  results  obtained  for  poration  in  normal  and  cancer  cells  due  to 
different  pulses  at  both  the  inner  and  outer  membranes  are  given  in  Table  I.  Basically, 
the  following  trends  are  evident: 

(i)  It  appears  on  the  basis  of  the  present  results  that  the  inner  membrane  of  a  normal, 
healthy  cell  is  difficult  to  porate. 

(ii)  With  increasing  electric  field,  the  inner  organelle  membrane  can  be  porated.  This 
should  lead  to  possible  apoptosis  or  cell  death  due  to  such  ultra-short,  high- 
intensity  pulses. 

(iii)  Longer  duration  pulses  are  more  likely  to  do  more  damage.  For  example,  at  the 
100  kV/cm  field,  the  50  ns  pulse  is  predicted  to  porate  both  the  inner  and  outer 
membranes  of  the  cancer  B-cell,  while  a  10  ns  duration  would  leave  the  inner 
membrane  intact. 


Table  I 


Normal  B-cell 

Outer 

membrane 

Inner 

Membrane 

E:  15kv/cm, 

No 

No 

E:  45kv/cm, 

Yes 

No 

E:  lOOkv/cm, 

Yes 

No 

Cancer  R-cell 

E:  45kv/cm, 

Yes 

No 

E:  lOOkv/cm. 

Yes 

No 

E:  45kv/cm, 

Yes 

No 

E:  lOOkv/cm, 

Yes 

Yes 

E:  225kv/cm, 

1 .5ns+10ns+1 .5ns 

Yes 

Yes 

The  Dynamic  Pore  Model: 

Coupling  the  above  electrical  model  with  the  Smoluchowski  equation  (SE)  for  the 
growth  and  decay  of  pores  leads  to  a  dynamic  simulation  of  the  cell  to  an  electrical  pulse. 
The  concept  of  using  a  diffusive  motion  across  an  energy  landscape  was  originally 
proposed  by  Smoluchowski  [8],  reviewed  by  Chandrasekhar  [9],  and  subsequently  used 
by  Kramers  to  develop  the  theory  of  activated  transition  rates  in  the  limit  of  large  friction. 
As  a  consequence  of  the  fluctuating-dissipative  nature  of  the  bio-system,  the  evolution  of 
the  mechanical  variables  (the  pore  density  in  this  case)  can  be  expected  to  follow  a 
diffusive  Brownian  motion  across  the  energy  landscape.  The  continuum  Smoluchowski 
theory  yields  the  following  for  the  pore  density  distribution  function  n(r,  t),  with  r  being 
the  pore  radius  [11,12]: 


dnjrj)  D  d(«(r,t)3E(r)/0r)  d2n(r,t ) 

dt  kBT  dr  dr2 


(5) 


where  S(r)  is  the  source  (or  pore  formation)  term,  while  D  is  a  pore  diffusion  constant, 
and  E(r)  is  the  pore-formation  energy.  Details  of  the  SE,  its  application  and  the 
dynamical  aspects  of  electroporation  have  been  addressed  in  a  series  of  recent  papers 
[1,5,6]  by  our  group,  and  will  not  be  reviewed  here.  The  pore  area  Ap  is  governed  by  the 
pore  density  distribution  function  n(r,t),  which  in  turn,  obeys  the  Smoluchowski  equation. 
It  turns  out  that  the  dynamic  pore  area  Ap  is  given  by: 


Ap(r,t)  ~  \[  1 2nrn(r\t)dr]  ,  (6) 

where  A0  is  the  equilibrium  area  under  zero  bias.  As  pores  form  in  response  to  an 
externally  applied  electric  pulse,  the  value  of  AP  increases,  and  the  electrical  conductivity 
through  the  pore  becomes  time-varying.  The  effective  membrane  conductivity 
aeff  _mem(0 can  be  calculated  by: 


eff  _  mem 


(0  = 


AJr,t) 


+  0- 


AP(r,t ) 


(7) 


where  <7  is  the  conductivity  of  the  medium  outside  the  membrane  and  rr  is  the 

conductivity  of  the  unporated  membrane.  Thus,  the  SE  yields  both  the  pore  distribution 
and  dynamic  pore  area,  which  in  turn  affects  the  electrical  parameters.  This  is  coupled  to 
the  distributed  electrical  analysis  for  self-consistency  through  dynamic  variations  of  the 
resistance  R  and  capacitance  C. 


(B)  Energy  Landscape  Calculations  for  Pulse-Width  Dependence: 

The  above  macroscopic,  distributed  electrical  model  was  applied  to  probe  and 
analyze  the  pulse-width  dependence  of  cellular  response.  Experimental  data  indicates 
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that  for  a  given  fixed  pulse  energy,  cell  apoptosis  is  more  pronounced  for  longer  pulses 
[13].  Some  of  the  experimental  results  are  shown  in  Fig.  4.  The  data  shows  caspase 
activity  (an  apoptosis  marker)  in  Jurkat  cells  subjected  to  a  single  pulse  of  variable 
duration,  but  fixed  total  energy  of  1.7  J/cc.  Longer  duration  pulses  are  again  seen  to 
produce  stronger  apoptotic  responses. 

We  have  attempted  to  provide  a  simple,  but  physical  model  for  cell  irreversibility 
and  apoptosis  in  response  to  an  ultrashort  (nanosecond),  high-intensity  electric  pulse.  We 
hypothesize  that  apoptotic  behavior  may  be  the  result  of  irreversible  conformational 
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Fig.  4.  Data  showing  caspase  activity  in  Jurkat  cells  subjected  to  a  single  pulse  of 
variable  duration,  but  fixed  total  energy  of  1.7  J/cc.  Longer  duration  pulses  are  again 
seen  to  produce  stronger  apoptotic  response. 

changes  at  the  inner  mitochondrial  membrane.  Such  changes  are  likely  to  be  driven  by 
the  high  electric  fields  that  arise  from  strong  transient  increases  in  transmembrane 
potentials.  Our  approach  was  based  on  an  energy  landscape  model,  which  is  used  to 
determine  the  temporal  evolution  of  the  configurational  probability  distribution  function 
[5].  Such  an  energy  landscape  picture  has  been  used  in  the  description  of  protein  changes 
[14],  as  well  as  the  folding/unfolding  of  small  globular  proteins  [15].  A  comprehensive 
and  accurate  treatment  of  this  problem  is  extremely  difficult  and  challenging  given  that 
apoptosis  can  have  several  pathways,  the  sequential  details  of  the  bio-physics  are  not  well 
known,  and  neither  are  the  magnitudes  of  the  internal  energies  and  configurations.  Given 


the  various  difficulties,  the  present  model  merely  attempted  to  provide  a  possible  physical 
picture  dependent  only  on  the  electrical  trigger,  and  captured  the  experimentally  observed 
pulse-width  dependence. 

Our  model  was  probabilistic  through  the  use  of  the  Smoluchowski  equation  and,  it 
coupled  the  distributed  electrical  model  for  current  flow  to  provide  the  time-dependent 
trans-membrane  voltages.  The  results  are  shown  in  Fig.  5,  and  agree  with  the  observed 
experimental  data  very  well.  The  probability  for  occupancy  of  a  secondary  energy 
minima  within  a  complex  energy  landscape  is  shown.  The  initial  state  of  the  system  has 
complete  and  total  occupancy  of  the  lowest  energy  valley.  With  energy  supplied  from 
the  external  electric-pulse,  the  occupancy  shifts  to  higher  energy,  and  some  states  in  the 
secondary  valley  can  be  occupied.  The  barrier  between  the  two  adjacent  valleys 
precludes  complete  recovery,  and  hence,  an  inherent  irreversibility  within  the  system 
results. 


Figure  5a.  Simulation  results  of  the  occupancy  distribution  within  a  two-valley,  energy 
landscape  model,  (i)  Snapshots  of  the  equilibrium  function  and  at  various  times 
following  10  ns,  60  ns  and  300  ns  pulses. 

Based  on  the  model  developed,  we  predict  a  number  of  interesting  features. 

(i)  First,  cell  irreversibility  at  a  fixed  input  energy  will  depend  on  the  pulse  width, 
and  is  likely  to  have  an  optimal  range.  Pulses  that  are  very  short  would  not  have  a 
significant  effect  due  to  insufficient  durations  for  electric  field  driven 
conformational  changes.  On  the  other  hand,  much  slower  and  longer  pulses 
would  also  be  rendered  ineffective  due  to  their  inability  to  penetrate  the 
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membranes  of  the  inner  organnelles  and  develop  a  significant  transmembrane 
voltage. 


Normalized  q/q1 

Figure  5b.  Simulation  results  of  the  occupancy  distribution  within  a  two-valley,  energy 
landscape  model.  Plots  emphasizing  higher  configuration  p(q)  values  within  the 
secondary  valley. 


(ii)  The  electrical  effects  are  likely  to  be  cumulative  and  to  exhibit  a  memory  effect, 
provided  successive  pulses  were  to  be  applied  with  delay  times  less  than  that 
required  for  complete  recovery  of  the  p(q)  profiles.  Such  a  memory  effect  has 
indeed  been  noted  in  the  context  of  the  ultrashort,  high-intensity  pulsing. 

(iii)  The  inherent  probabilistic  basis  of  this  model  implies  that  complete  and  total 
killing  of  cell  populations  can  never  be  achieved  by  a  single  shot  exposure.  There 
will  always  be  a  finite  probability  distribution  below  the  threshold  because  of  the 
diffusive  motion  in  energy  space.  Only  at  very  low  temperatures  might  one 
expect  cohesive  and  “ballistic”  motion.  This  prediction  of  incomplete  killing  is  in 
keeping  with  experimental  observations. 

(iv)  Furthermore,  this  model  predicts  that  the  use  of  very  high  electric  field  intensities 
may  not  necessarily  have  a  significantly  larger  impact.  This  naturally  follows 
since  not  only  is  the  pulse  duration  an  important  factor,  but  also  because  the 
transmembrane  potential  would  not  scale  with  the  applied  field.  Electroporation 
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and  dynamic  changes  in  conductivity  would  work  to  offset  the  impact  of  larger 
external  fields. 

(v)  Here  the  irreversibility  mechanism  is  taken  to  arise  from  fixed  charge  and  dipole 
movements  brought  about  within  the  membrane  by  the  electrostatic  driving  force 
produced  by  the  large  transient  transmembrane  potentials.  Consequently, 
excitation  processes  that  cannot  generate  large  electric  fields  at  membranes,  or  if 
their  time  durations  are  small,  might  not  be  very  effective  in  electrically  triggering 
apoptotic  behavior  in  cells.  For  example,  under  “contactless”  conditions,  or  for 
microwave  and  radio-frequency  excitation,  cellular  apoptosis  may  be  expected  to 
be  fairly  negligible. 

(C)  Implementation  of  Molecular  Dynamics  Simulations  for  Membrane  Response: 

Microscopic  analysis  based  on  the  Molecular  Dynamics  atomistic  approach  was 
developed  and  applied  to  probe  the  membrane  response.  The  microscopic  MD  simulation 
approach  is  superior  to  continuum  approaches.  Advantages  of  this  approach  are:  (i) 
Inclusion  of  the  collective,  many-body  interaction  potentials  at  the  nano-scale  level,  (ii) 
Dynamical  screening,  (iii)  Avoidance  of  the  “mean-field”  approximations,  (iv)  A  natural 
inclusion  of  noise  and  statistical  fluctuations,  (v)  Self-consistent  and  dynamical  transport 
calculations  without  arbitrary  fitting  parameters,  and  (vi)  Easy  incorporation  of  arbitrary 
defects  and  non-uniformities,  as  well  as  complex  geometries. 

In  our  studies,  the  GROMACS  package  [16,17]  was  used  for  the  MD  simulations 
of  field-induced  membrane  effects.  It  was  applied  to  dipalmitoyl-phosphatidyl-choline 
(DPPC)  bilayers.  Details  of  the  dode-cylpho-sphocholine  (DPC)  force  field  for 
membrane  molecular  motion  can  be  found  elsewhere  [18].  The  system  was  coupled  to  an 
isotropic  pressure  bath  of  1  atmosphere,  and  a  heat  bath  of  300  K  [19].  Velocities  of  the 
water  and  membrane  molecules  were  scaled  with  a  coupling  constant  of  0.1  ps.  The 
algorithms  for  pressure  and  temperature  control  were  those  discussed  by  Allen  and 
Tildesly  [20].  A  2  fs  time  step  was  used  with  constrains  on  the  bond  lengths  within  the 
lipids  and  on  the  water  geometry.  A  group  based  twin  cutoff  scheme  was  employed  for 
the  non-bonded  interactions,  with  Cutoff  of  1.0  nm  for  Lennard-Jones,  and  RC[S toff  of  1.5 
nm  for  electrostatic  interactions.  The  SPC  water  model  was  used  with  a  simple  point 
charge  of  -0.82e  at  the  oxygen  atom,  and  charges  of  0.4  le  on  each  of  the  hydrogen 
atoms.  The  O-H  bond  length  was  0.1  nm  with  a  109"  28'  angle  for  the  H-O-H  bond. 
Both  the  bond  length  and  bond  angle  for  water  are  constrained  by  the  SHAKE  method 
[21].  All  ions  were  treated  as  single  Lennard-Jones  centers,  with  the  parameters  for 
energy  (ejj)  and  distance  (oy)  taken  from  published  literature  [22,16,17].  Periodic 
boundaries  were  applied,  and  the  simulations  carried  out  under  constant  pressure 
conditions.  This  method  [23],  allows  the  simulation  box  size  to  change  so  that  the 
internal  virial  matches  the  external  pressure. 

The  cell  membrane  was  taken  to  comprise  of  Dioleoyl-phosphatidyl-choline 
(DPPC)  lipid  molecules.  These  are  complex  chain  containing  a  large  number  of  atoms, 
with  a  prominent  “head”  and  two  “tails”  as  shown  in  Fig.  6a.  In  order  to  ease  the 
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computational  burden  and  facilitate  calculations  in  a  reasonable  amount  of  time,  the 
complete  DPPC  molecule  was  reduced  to  a  coarse  grained  model  as  shown  in  Fig.  6b. 
Thus,  the  choline,  phosphate,  glycerol,  and  methyl  groups  within  DPPC  were  each 
aggregated  to  just  one  mega-molecule  of  large  mass.  The  mega-molecules  were  treated 
as  single  Lennard-Jones  centers,  with  the  energy  and  distance  parameters  for  each  sub¬ 
group  taken  from  the  literature  [16,17,22].  Such  coarse-grained  models  have  been  used 
recently  [24]  to  reduce  the  computational  burden  associated  with  an  all-atom  approach 
[25], 


Fig. 6.  The  dipalmitoyl  -phosphatidyl-choline  (DPPC)  lipid  molecule,  (a)  All  atom 
representation,  and  (b)  the  coarse-grained  model  used. 

For  the  actual  simulation,  a  total  of  128  DPPC  chains  (arranged  on  a  uniform  8  x 
8  square  grid)  and  3655  water  molecules  were  used,  yielding  a  ratio  of  28.56: 1  between 
water  and  the  lipids.  The  total  simulation  volume  was  taken  to  a  6.8  nm  x  6.8  nm  x  6.8 
nm  cubic  box.  The  lipid  bilayer  was  taken  to  have  an  initial  thickness  of  4.5  nm.  The 
starting  configuration  was  created  by  assigning  the  DPPC  molecules  at  each  of  the  64 
locations  within  the  square  grid.  The  water  molecules  were  randomly  placed  within  the 
volume  lying  at  the  top  and  on  the  bottom  of  the  DPPC  layer.  Due  to  our  use  of  a 
constant  pressure  coupling,  the  initial  size  is  not  critical  since  the  system  continually 
adjusts  to  match  the  external  pressure.  In  all  cases,  before  starting  the  field-dependent 
calculations  of  interest,  this  system  was  allowed  to  reach  its  self-consistent  energy 
minimized  state  by  running  the  MD  simulations  at  zero  electric-field  for  100  ps 

Verification  and  validation  of  the  MD  code  was  first  carried  out,  and  the  results 
compared  with  available  reports.  Towards  this  goal,  simulations  at  zero  electric  field 
were  first  carried  out  in  keeping  with  a  recent  publication  by  Lopez  et  al.  [24].  Under 
this  zero-field  conditions,  one  does  not  expect  pores  to  form,  or  for  water  molecules  to 
pass  through  the  membrane.  Our  results  were  in  agreement  with  expectations  and 
demonstrated  this  qualitative  trend.  This  is  clearly  seen  from  the  spatial  electron  density 
distribution  result  shown  in  Fig.  7.  The  absence  of  electron  density  over  the  central 
portion  of  the  membrane  is  indicative  of  the  inability  of  water  molecules  to  pass  through 
the  DPPC  membrane. 
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Electron  Density 


Figure  7.  Molecular  Dynamics  simulation  results  showing  the  electron  density  distribution 
at  zero  electric  Field.  Water  molecules  do  not  penetrate  the  central  membrane  region. 


Next,  the  results  of  applying  a  realistic  high-intensity,  ultra-short  electric  pulse  are 
presented  and  discussed.  Snapshots  showing  the  process  of  membrane  configurational 
change  are  given  in  Fig.  8.  Figure  8(a)  shows  the  membrane  to  be  well  aligned 
(designated  as  the  t  =  0  ns  instant)  right  after  the  energy  minimization  was  performed. 
After  about  4.8  ns,  some  structural  rearrangement  (i.e.  defect  formation)  begins  to  form  at 
the  anode  side  (i.e.,  the  membrane  bottom)  as  shown  in  Fig.  8(b).  A  pore  is  then  seen  to 
form  fairly  quickly  after  the  creation  of  the  initial  defect.  In  fact,  a  complete  nano-pore 
formation  is  predicted  within  about  0.5  ns  after  the  initial  structural  change.  Obviously, 
most  of  the  time  for  the  electroporation  process  is  taken  up  by  the  initial  defect  formation. 
A  pore,  formed  at  5.3  ns,  is  shown  in  Fig.  8(c). 

This  result  demonstrates  several  important  points.  First,  there  is  a  very  good  match 
between  the  MD  result  of  Fig.  8,  predictions  of  our  distributed  current  model,  and  actual 
experimental  data.  In  all  cases,  a  nano-pore  is  predicted  to  occur  within  about  5-6  ns  at 
these  high  field  intensities  applied  to  the  cells.  Since  the  pulse  duration  is  typically 
longer  than  this  time,  ion  transport  (especially  for  ions  with  smaller  radii)  would  be  very 
likely.  Second,  the  initial  structural  rearrangement  and  dipole  re-orientation  is  a  critical 
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step  in  the  electroporation  process.  Once  such  an  initial  break-through  is  achieved,  the 
poration  process  proceeds  relatively  quickly.  Finally,  the  poration  process  has  a  polarity 


Fig.  8.  MD  snapshots  of  the  electroporation  process  in  a  DPPC  membrane,  (a)  Initial 
configuration  at  t  =  0  ns;  (b)  Structural  re-arrangement  begins  at  the  membrane  bottom 
(anode  side)  at  about  t  =  4.8  ns;  (c)  A  pore  predicted  to  form  at  t  =  5.3  ns. 

dependence,  and  begins  on  the  anodic  side  of  a  membrane.  The  related  physics  can 
easily  be  understood  by  considering  the  dipolar  configuration  within  the  DPPC 
membrane.  Fig.  9  shows  a  simple  schematic  of  the  membrane  lipids  with  their  dipoles 
located  at  the  head-groups.  For  each  DPPC  chain,  the  head  group  contains  a  dipole  with 
positive  charge  on  Choline  and  the  negative  charge  centered  on  the  Phosphate  group. 
Initially,  with  no  electric  field  present,  the  dipoles  are  in  random  thermal  motion  with  the 
positive  charges  residing  on  the  outermost  portions  of  the  lipid.  Electrical  field  induced 
defects  are  initiated  by  the  movement  of  dipoles  on  the  surface  of  the  membrane.  Defects 
start  to  form  on  the  anode  side  of  the  membrane  because  positively  charged  molecules 
(e.g.,  Choline)  on  this  side  are  forced  to  swing  around  (i.e.,  re-orient  in  the  presence  of  a 
strong  external  electric  field)  and  enter  the  membrane.  However,  this  same  applied 
electrical  field,  when  acting  on  the  dipoles  located  at  the  cathodic  membrane  surface 
merely  works  to  stretch  the  dipoles  without  any  molecular  movement  into  the  membrane 
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volume.  Thus,  Fig.  9(a)  shows  randomly  distributed  dipoles  at  the  head  groups  on  either 
side  of  the  membrane  as  the  initial  configuration.  Fig.  9(b)  shows  the  alignment  of 


Fig.  9.  Schematic  of  the  pore  initiation  process  at  the  DPPC  membrane.  The  arrows 
present  dipoles  on  the  head  groups  of  lipid  molecules,  (a)  Initial  configuration  without  an 
externally  applied  field,  and  (b)  the  configuration  and  expected  movement  of  the  dipole 
head  groups  upon  the  application  of  an  external  electric  field. 


dipolar  head  groups  on  the  anodic  side  gradually  deviating  from  the  normal  (equilibrium) 
orientation,  and  a  defect  starting  to  form. 


If  PS  externalization  is  a  pore-driven  event,  then  this  process  should  also  start  at 
the  anode-side.  Experimental  observations  do  indicate  just  such  an  anode-side  preference. 
The  electrostatics  provides  an  additional  rationale  for  an  anode-side  event.  Since  the  PS 
molecules  residing  on  the  inner  leaflet  are  negatively  charged,  the  externally  applied 
electric  field  will  tend  to  push  out  PS  at  the  anode  side,  while  pulling  it  inwards  on  the 
cathodic  side.  Coupling  this  with  an  anode-side  pore  formation  event,  leads  to  the 
collective  effect  of  preferential  PS  externalization  at  the  anode  end.  Fig.  10  is  a  simple 


Fig.  10.  Schematic  of  the  pore  initiation  process  within  a  DPPC-DPPS  membrane.  The 
arrows  represent  dipoles  on  the  lipid  head  groups,  and  the  shaded  circle  denotes  the 
negative  DPPS  head  group,  (a)  Initial  configuration  without  an  externally  applied  field; 
(b)  configuration  of  the  dipole  head  groups  with  an  external  field;  (c)  Nano-pore 
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formation  with  a  DPPS  molecule  close  to  the  pore;  and  (d)  DPPS  molecular  movement  to 
the  outer  leaflet  of  the  membrane. 

schematic  demonstrating  the  dynamics.  Initially  (Fig.  10a),  the  dipoles  at  the  DPPC 
lipids  (shown  as  arrows)  are  randomly  located,  with  a  negative  PS  on  the  inner  leaflet. 
With  application  of  the  electric  field,  the  dipoles  on  the  outer  layer  (anode-side)  re-orient 
leading  to  defect  initiation  at  the  outer  membrane  surface.  Eventually  a  pore  forms  (Fig. 
10c),  and  the  negatively  charged  PS  begins  to  drift  and  diffuse  towards  the  exterior 
surface.  Results  of  MD  simulations,  shown  in  Figs.  1  la- 1  Id  demonstrate  the  above  more 
clearly.  An  initial  snapshot  of  the  membrane  system  is  shown  in  Fig.  1 1(a),  with  the  PS 
molecule  located  on  one  side  of  the  membrane,  opposite  the  anode.  Fig.  1 1(b)  shows  a 
pore  starting  to  form  with  some  translocation  of  the  PS  chain  at  3.2  ns.  Due  to  a  large 
electrostatic  force  on  the  PS  head  group  (in  the  range  of  _10-10N),  the  chain  is  dragged 
half  way  to  the  anode  side  of  the  membrane  along  the  wall  of  the  nanopore  at  about  _3-34 
ns  as  shown  in  Fig.  1 1(c).  Finally,  in  Fig.  1 1(d),  the  DPPS  chain  is  on  the  other  leaflet  of 
the  membrane  at  _3.61  ns.  The  MD  simulation  thus  validates  the  pore-facilitated,  field- 
assisted  mechanism  of  PS  externalization. 


Fig.  1 1.  Snapshots  of  the  PS  externalization  process  at  a  DPPC  membrane,  (a)  Initial 
configuration  at  _  0  ns;  (b)  A  pore  forms  and  DPPS  starts  to  go  along  the  pore  wall  at 
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3.2  ns;  (c)  DPPS  half  way  to  externalization  at  _3.34ns;  and  (d)  DPPS  on  the  other  leaflet 
of  the  membrane  at  _3.61  ns. 

Objective  II.  Time  Domain  Dielectric  Spectroscopy 

The  second  objective  of  the  Bio-Inspired  Concepts  grant  focused  on  setting  up 
and  developing  a  state-of-the-art  diagnostic  system  for  monitoring  the  electrical 
properties  of  biological  cells.  It  was  based  on  utilizing  Time  Domain  Dielectric 
Spectroscopy  (TDDS)  that  allows  for  the  measurement  of  frequency-dependent  complex 
permittivity  and  conductance  parameters  of  biological  cells  and  tissues. 

Accomplishments  Under  Objective  II 

Pulsed  electric  fields  (PEFs)  above  a  certain  voltage  threshold  cause  electroporation  for 
microsecond  pulses  and  intracellular  effects  for  sub-microsecond  pulses.  Models 
describing  these  effects  often  depend  on  the  electrical  properties  of  the  cell,  which  are 
altered  by  the  PEF.  We  have  used  time  domain  dielectric  spectroscopy  (TDDS)  to 
provide  data  for  these  models  and  to  measure  effects  on  cell  suspension  conductivity.  In 
particular,  single  50  ps  and  10  ns  pulses  were  applied  to  HL-60  cells,  with  the  voltages 
chosen  so  the  pulses  have  approximately  the  same  energy.  For  1.1  kV/cm,  50  ps  pulses, 
the  conductivity  rose  within  a  minute  after  the  pulse  and  dropped  dramatically 
approximately  forty  minutes  after  the  pulse.  For  78  kV/cm,  10  ns  pulses,  a  brief  delay 
prior  to  the  conductivity  rise  was  observed  and  the  same  drop  in  conductivity  after 
approximately  forty  minutes  was  noted.  For  both  pulse  durations,  higher  voltages 
frequently  led  to  membrane  poration  followed  by  a  gradual  recovery  approximately  thirty 
to  forty  minutes  after  the  pulse.  For  2  kV/cm,  50  ps  pulses,  significantly  larger  scatter  in 
Trypan  Blue  uptake  measurements  were  observed  due  to  stronger  effects  on  the  cell 
membrane.  By  using  a  cell  model,  it  was  shown  that  a  50  ps  pulse  caused  a  much  larger 
rise  in  membrane  conductivity  than  a  10  ns  pulse  of  the  same  energy. 

Figure  12  shows  a  block  diagram  of  the  Time  Domain  Spectrometer  (IDC  Expertise, 
United  Kingdom)  based  on  a  design  by  Feldman’s  group  [26,  27].  An  external  pulse 
generator  provides  a  200  mV,  5  ps  pulse  with  a  40-50  ps  rise-time.  The  cell  suspension 
can  be  placed  in  either  a  coaxial  (high  frequency)  sample  holder  or  a  parallel  plate  (low 
frequency)  sample  holder  [27].  The  coaxial  sample  holder  was  chosen  because  it  allowed 
for  the  use  of  fewer  cells  while  still  allowing  measurements  in  the  frequency  range  of 
interest.  Due  to  the  strong  temperature  dependence  of  the  dielectric  parameters,  a 
thermostat  (Julabo,  US)  was  used  to  maintain  sample  cell  at  25  C.  After  suitable  signal 
processing  [26,  27],  a  Fourier  Transform  was  used  to  obtain  the  real  and  imaginary 
components  of  permittivity  for  the  cell  suspension.  From  this  data,  the  permittivity  and 
conductivity  of  the  cell  membrane,  cytoplasm,  nuclear  envelope,  and  nucleoplasm  could 
all  be  obtained  for  HL-60  using  the  two-shell  model  [26,  28]. 
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Figure  12.  Block  diagram  of  the  TDDS  system. 

The  results  for  the  time-dependent  permittivity  (em)  obtained  are  shown  in  Fig.  13 
after  a  typical  140  kV/cm,  10  ns  pulse  and  after  a  typical  2  kV/cm,  50  ps  pulse.  For  the 
140  kV/cm  pulse,  €m  rises  from  an  initial  control  value  of  5.70  to  8.02  within  ninety 
seconds  from  the  pulse,  and  then  decreases  to  a  5.345  value  fourteen  minutes  after  the 
pulse.  For  the  2  kV/cm  pulse,  em  decreased  from  3.04  to  2.35  thirty-five  seconds  after  the 
pulse,  and  then  rose  to  3.32  about  nineteen  minutes  after  the  pulse. 


Fig.  1 3.  Membrane  permittivity  (em)  of  an  HL-60  cell  as  a  function  of  time  after  a  2 
kV/cm,  50  ps  pulse  and  a  140  kV/cm,  10  ns  pulse. 
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Figure  14  shows  the  conductivity  (om)  of  an  HL-60  cell  membrane  as  a  function  of  time 
after  a  140  kV/cm,  10  ns  pulse  and  after  a  2  kV/cm,  50  ps  pulse.  For  the  140  kV/cm 
pulse,  am  rose  from  7.5  x  10'6  S/m  to  1.2  x  10'^  S/m  about  ninety  seconds  after  the  pulse. 
It  then  was  seen  to  decrease  to  6.7  x  10"6  S/m  about  fourteen  minutes  after  the  pulse.  For 
2  kV/cm  pulse,  am  rose  from  3.6  x  10  6  S/m  to  1.08  x  10'5  S/m  thirty-five  seconds  after 
the  pulse  and  9  x  10‘6  S/m  about  nineteen  minutes  after  the  pulse. 


Time  after  pulse  (min) 

Fig.  14.  Membrane  conductivity  (am)  of  an  HL-60  cell  membrane  as  a  function  of  time 
after  a  140  kV/cm,  10  ns  pulse  and  a  2  kV/cm,  50  ps  pulse. 

The  next  step  in  this  analysis  is  to  utilize  the  full  two-shell  model  to  determine  the 
electrical  properties  of  each  component  of  the  cell  [28].  A  complete  and  rigorous 
analysis  has  so  far  been  difficult  for  a  variety  of  reasons.  One  problem  is  that  the  two- 
shell  model  is  very  sensitive  to  the  geometric  parameters  of  the  cell  [28],  which  can  be 
greatly  altered  by  the  PEF.  Moreover,  electroporation  and  intracellular  manipulation  are 
inherently  anisotropic. 

Objective  III.  Cellular  Measurements  of  High-Field  Pulse  Effects 

A  final  portion  of  the  Bio-Inspired  Concepts  grant  focused  on  experimental 
measurements  of  cell  response  to  pulsed  electric  fields. 

Accomplishments  Under  Objective  III 

The  effects  of  nanosecond  (10-300ns)  pulsed  electric  fields  (nsPEF)  on 
mammalian  cell  structure  and  function  were  probed.  As  the  pulse  durations  decreased, 
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effects  on  the  plasma  membrane  (PM)  correspondingly  decreased  and  effects  on 
intracellular  signal  transduction  mechanisms  increased.  When  nsPEF-induced  PM 
electroporation  effects  did  occur,  they  were  distinct  from  classical  PM  electroporation 
effects,  suggesting  unique,  nsPEF-induced  PM  modulations.  NsPEF  with  durations  and 
electric  field  intensities  were  seen  to  induce  apoptosis  in  mammalian  cells  with  a  well- 
characterized  phenotype  typified  by  extemalization  of  phosphatidylserine  on  the  outer 
PM  and  activation  of  caspase  proteases. 

Some  of  the  relevant  experimental  details  are  summarized  below. 

Cell  culture:  HL-60  and  Jurkat  cells  were  obtained  from  and  cultured  as  recommended  by 
American  Type  Culture  Collection  as  described  in  detail  elsewhere  [13,  29,  30].  Cells 
were  removed  from  log  phase  growth  and  suspended  in  Hanks  Balanced  Salt  Solution 
(HBSS)  or  Tyrodes  buffer. 

Administration  of  nsPEF:  Cell  suspensions  (7.7  x  106  cells/ml;  106  cells/1 30pl)  were 
loaded  into  the  BioRad  gene  Pulser®  cuvettes  (Bio-Rad  laboratories;  Hercules,  CA)  prior 
to  nsPEF  pulsing.  The  pulse  was  delivered  to  the  load,  via  either  a  0.1,  0.2,  or  0.4  cm 
cuvette  containing  cell  suspension,  by  means  of  a  cable  pulse  generator.  Briefly,  the 
generator  consists  of  a  pulse-forming  network  (PFN),  -  five  50  cables  in  parallel  and  a 
spark  gap  in  atmospheric  air  as  a  nanosecond  closing  switch.  Post  pulse  the  cell 
suspension  was  removed  from  the  pulsing  cuvette  and  assayed. 

Determination  of  plasma  membrane  integrity:  HL-60  or  Jurkat  cells  were  exposed  to 
nsPEF  in  the  presence  of  ethidium  homodimer- 1  or  propidium  iodide  (Molecular  Probes, 
Eugene,  OR)  and  analyzed  by  flow  cytometry.  Values  are  expressed  as  the  geometric 
mean  fluorescence  intensity  and/or  as  the  percentage  of  cells  that  exhibited  fluorescence. 

Effect  of  nsPEF  on  apoptosis  induction:  Apoptosis  was  determined  using  flow  cytometry 
and  well  defined  apoptosis  characteristics  including  maintenance  of  membrane  integrity 
as  determined  by  the  exclusion  of  ethidium  homodimer,  the  presence  of  annexin-V 
binding  as  a  marker  for  phosphatidylserine  extemalization  on  the  plasma  membrane  of 
intact  cells,  and  the  presence  of  caspase  activation  using  two  different  assays. 

Flow  cytometry:  Flow  cytometry  was  carried  out  as  previously  described.  Briefly,  ~10(' 
cells  were  exposed  to  nsPEF  in  electroporation  cuvettes.  Cell  suspensions  (1.5  x  104  /ml) 
were  then  analyzed  using  a  Becton-Dickinson  FACS  calibur  flow  cytometer. 
Fluorescence  emission  was  collected  with  interference  bandpass  filter  sets  for  (FL1-H) 
and  (FL2-H).  Data  was  stored  and  analyzed  by  CellQuest  software. 


Based  on  an  electric  model  of  biological  cells,  it  was  predicted  that  as  the  pulse 
duration  decreases  effects  on  the  plasma  membrane  would  decrease.  To  test  this 
hypothesis  directly,  Jurkat  or  HL-60  cells  were  exposed  to  pulsed  fields  with  durations  of 
10  or  60  ns  and  the  number  and  amplitude  of  the  pulses  were  varied  as  shown  in  figure 
15.  The  cells  were  treated  with  the  pulses  in  the  presence  of  ethidium  homodimer,  a 
fluorescent,  membrane  impermeable  DNA  stain  that  is  used  as  a  marker  for  plasma 
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membrane  integrity,  and  then  evaluated  for  fluorescence  by  flow  cytometry.  Panels  A 
and  B  show  results  from  Jurkat  cells  that  were  exposed  to  single  pulses  of  10  and  60  ns, 
respectively,  with  increasing  electric  fields.  The  level  of  fluorescence  is  a  direct  indicator 
of  levels  of  effects  on  plasma  membrane  electroporation.  For  any  given  electric  field, 
effects  on  plasma  membrane  integrity  were  less  for  the  shorter  pulses.  The  thresholds  for 
breaching  plasma  membrane  integrity  were  between  150  and  300  kV  cm"1  for  10  ns 
pulses  and  between  26  and  40  kV  cm"1  for  60  ns  pulses.  Thus,  higher  electric  fields  were 
required  to  affect  the  plasma  membrane  for  10  ns  versus  60  ns  pulses.  Similar  effects 
were  observed  for  HL-60  cells  (data  not  shown).  Panels  C  and  D  show  results  from  HL- 
60  cells  for  10  and  60  ns  pulses  with  the  electric  field  adjusted  to  provide  a  similar  energy 
density  with  increasing  pulse  numbers.  Effects  on  the  plasma  membrane  increased  with 
increasing  pulse  number  for  both  10  and  60  ns  pulses,  but  effects  at  10  ns  were  less  than 
the  effects  at  60  ns  for  any  given  pulse  number. 
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Figure  15.  Data  showing  that  as  the  pulse  number,  durations  and/or  the  electric  field 
intensities  decrease,  effects  on  the  plasma  membrane  (PM)  also  decrease. 
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In  order  to  compare  the  effects  of  nsPEF  and  classical  electroporation  on  plasma 
membrane  integrity,  Jurkat  cells  were  exposed  to  a  single  10  ns  pulses  at  150  kV  cm'1 
shown  in  figure  16  (panel  A)  or,  a  single  6  ms  pulse  at  0.5  kV  cm'1  in  the  presence  of 
ethidium  homodimer  and  analyzed  by  flow  cytometry.  The  effects  of  the  10  ns  pulse  on 
plasma  membrane  integrity  were  minimal.  In  contrast,  the  classical  electroporation 
pulses,  which  provided  an  energy  density  50%  lower  than  the  10  ns  pulse,  resulted  in 
considerable  increases  in  ethidium  homodimer  fluorescence,  indicating  significant 
disruption  of  the  plasma  membrane.  When  ethidium  homodimer  was  added  10  or  45  min 
after  exposure  to  each  of  the  pulsed  electric  fields,  no  increase  in  fluoresce  was  observed 
after  the  10  ns  pulse,  but  significant  fluorescence  was  observed  after  the  classical 
electroporation  pulses  (data  not  shown).  These  results  indicate  that  nsPEF  have  minimal 
effects  on  the  plasma  membrane  that  are  short  lived,  while  classical  plasma  membrane 
electroporation  pulses  have  considerable  longer  lasting  effects  on  the  plasma  membrane 
and  are  thus  functionally  different. 
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Figure  16.  Classical  plasma  membrane  electroporation  and  nsPEF  conditions  have 
significantly  different  effects  on  plasma  membrane  integrity. 


As  illustrated  in  figure  17,  HL-60  cells  were  exposed  to  nsPEF  with  durations  of 
10  ns  (top  panels)  or  60  ns  (bottom  panels)  and  increasing  electric  fields  as  indicated. 
Cells  were  exposed  to  nsPEF  in  the  presence  of  ethidium  homodimer  as  a  marker  for 
plasma  membrane  integrity  (left  panels)  during  the  pulses  and  then  Annexin-V-FITC  was 
added  as  an  apoptosis  marker  for  PS  externalization  (right  panels).  Cells  were  analyzed 
for  fluorescence  by  flow  cytometry  10  min  post-pulse.  As  indicated  by  ethidium 
homodimer  fluorescence  (left  panels)  none  of  the  conditions  led  to  significant  losses  in 
plasma  membrane  integrity.  In  these  experiments,  no  significant  losses  of  plasma 
membrane  integrity  occurred;  only  3%  of  the  cells  exhibited  minimal  ethidium 
fluorescence  at  60  ns  and  60  kV  cm'1.  In  contrast,  Annexin-V-FITC  binding  was 
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observed  at  10  ns  with  electric  fields  of  greater  than  150  kV  cm'1  ,  and  at  60  ns  with 
electric  fields  of  greater  than  26  kV  cm'1.  Like  nsPEF-induced  effects  on  plasma 
membrane  integrity,  lesser  effects  were  observed  for  PS  extemalization  at  shorter  pulse 
durations,  and  higher  electric  fields  were  required  to  induce  PS  extemalization  as  the 
pulse  duration  decreased.  Thus,  it  was  possible  to  observe  nsPEF-induced  PS 
extemalization  under  conditions  that  did  not  include  loss  of  plasma  membrane  integrity. 
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Figure  17.  Phosphatidylserine  (PS)  extemalization  occurs  independent  of  loss  of  plasma 

membrane  integrity  in  response  to  nsPEF. 


SUMMARY 

A  comprehensive  modeling  and  experimental  effort  was  carried  out  to  develop  an 
understanding  of  cellular  bio-response  to  short  duration,  high-intensity  electric  fields.  The 
simulation  capability  included  macroscopic  models  for  determining  the  time-dependent 
spatially-variable  electric  potential  and  current  flows  at  single  cells.  This  provides  the 
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means  of  predicting  changes  in  the  transmembrane  voltages  at  the  plasma  and  other  intra¬ 
cellular  membranes.  This  voltage  change  is  thought  to  be  the  primary  driver  in  bringing 
about  change  (both  reversible  and  irreversible)  in  bio-systems.  The  macroscopic  model 
also  provides  for  predictions  of  heating  and  internal  temperature  rises. 

The  macroscopic  model  was  then  coupled  to  a  microscopic,  nano-simulator  to 
probe  the  sub-cellular  response  at  the  molecular  level.  Many  of  the  observed  details  such 
as  PS  extemalization,  the  time  scales  for  pore  formation,  and  their  probable  diameters 
were  predicted.  The  results  agreed  well  with  the  experimental  data,  and  provide  a  logical 
basis  for  the  lack  of  Propidium  Iodide  transport  during  cell  exposure  to  pulsed  voltages. 

For  accurate  predictions,  the  relevant  parameters  of  the  bio-system  need  to  be 
known  or  to  be  obtained.  This  was  done  within  our  group  by  developing  the  Time 
Domain  Dielectric  Spectroscopy  method.  The  system  is  operational,  and  yields  data  on 
the  conductivity  and  permittivities  of  cells  and  its  organelles.  This  is  a  useful  and 
important  development. 

Finally,  cellular  response  were  measured  based  on  a  variety  of  techniques, 
including  flow  cytometry,  optical  microscopy  and  imaging.  The  central  results  were:  (i) 
There  is  a  critical  electric  field  (or  transmembrane  potential)  below  which  no  cell  death  is 
observed.  Above  the  critical  field,  there  was  a  rapid  transition  from  100%  to  0%  cell 
survival,  (ii)  The  critical  membrane  voltage  for  irreversible  electroporation  decreases 
with  increasing  pulse  width  and  multiple  pulses  with  an  asymptotic  value  of  ~  1  volt  for 
long  pulse  lengths,  (iii)  Transmembrane  voltage  for  irreversible  electroporation  scaled 
with  constant  pulse  energy  for  pulse  lengths  less  than  5  (is.  (iv)  PS  extemalization 
leading  to  cell  death  through  macrophage  action  can  be  electrically  triggered.  Pore 
formation  can  be  reversible  at  the  short  nano-second  pulse  durations. 

This  study  would  be  indicative  of  the  following  conclusions:  (i)  It  may  be  energy 
efficient  to  use  short  pulses  for  cell  death.  The  elimination  of  heat  dissipation  in  the 
medium  would  contribute  to  this  efficiency,  (ii)  Selective  targeting  of  cells  for  apoptosis 
appears  to  be  possible  through  the  use  of  ultra-short  pulses,  (iii)  The  ability  to  extract 
energy  with  directed  motion  and  transport  from  low-amplitude  fluctuations  appears  to  be 
possible  in  cell  membranes  based  on  spatially  non-uniform  potentials. 
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